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Abstract 

We study high-p^ jets from QCD and from highly-boosted massive particles such as tops, W, Z and Higgs, and 
argue that infrared-safe observables can help reduce QCD backgrounds. Jets from QCD are characterized by different 
patterns of energy flow compared to the products of highly-boosted heavy particle decays, and we employ a variety 
of jet shapes, observables restricted to energy flow within a jet, to explore this difference. Results from Monte Carlo 
generators and arguments based on perturbation theory support the discriminating power of the shapes we refer to 
as planar flow and angularities. We emphasize that for massive jets, these and other observables can be analyzed 
perturbatively. 



Introduction. At the Large Hadron Collider (LHC), 
events with highly-boosted massive particles, tops W, 
Z and Higgs, h [2| may be the key ingredient for the dis- 
covery of physics beyond the standard model 0, 0, @] . In 
many decay channels, these particles would be identified 
as high-pr jets, and any such signal of definite mass must 
be disentangled from a large background of light-parton 
("QCD") jets with a continuous distribution. This back- 
ground far exceeds such signals, and relying solely on jet 
mass as a way to reject QCD background from signal 
would probably not suffice in most cases even using 
a narrow window for dijets in the search for massive par- 
ticles such as tops, produced in pairs. 

In this paper, we argue that for massive, high-py jets, 
infrared (IR) safe observables may offer an additional 
tool to distinguish heavy particle decays from QCD back- 
ground, perhaps even on an event-by-event basis. We will 
refer to inclusive observables dependent on energy flow 
within individual jets as jet shapes. A jet within a cone of 
radius 0.4, for example, may be identified from the energy 
recorded in roughly fifty 0.1x0.1 calorimeter towers. It 
thus contains a great deal of information. Perturbatively- 
calculable, infrared safe jet shapes combined, of course, 
with IR-safe jet finding algorithms [Tj, may enable us 
to access that information systematically, and to form a 
bridge between event generator output and direct theory 
predictions. 

Essentially, any observable that is a smooth functional 
of the distribution of energy flow among the cells defines 
an IR-safe jet shape [8]. The jet mass is one example, 
but a single jet may be analyzed according to a variety 
of shapes. In particular, the jet mass distribution has 



large corrections when the ratio of the jet mass to jet 
energy is small [9(, but can be computed at fixed order 
when the logarithm of that ratio has an absolute value of 
order unity. Once the jet mass is fixed at a high scale, a 
large class of other jet shapes become perturbatively cal- 
culable with nominally small corrections. Indeed, a jet 
whose mass exceeds O(100 GeV) becomes, from the point 
of view of perturbation theory, much like the final state in 
leptonic annihilation at a similar scale. At such energies, 
event shapes, which in the terminology of this paper are 
jet shapes extended over all particles, have been exten- 
sively studied in perturbation theory (Toj . In this study 
we explore the possibility that perturbative predictions 
for jet shapes differ between those jets that originate from 
the decay of heavy particles, and those which result from 
the showering of light quarks and gluons. Very interest- 
ing related studies have recently appeared in [4j, [ill, [l2j • 

Jet Shapes and Jet Substructure. We would like 
to identify jet shapes for which perturbative predictions 
differ significantly between QCD and other high-pr jets, 
focusing on relatively narrow windows in jet mass. In our 
companion paper [6( we have discussed how to calculate 
the jet mass distribution for the QCD background. We 
now extend this argument to the computation of other 
jet shape observables. 

We emphasize that, because the observables under con- 
sideration are IR-safe, we may calculate them as power 
series in a s , starting at first order for the QCD back- 
ground, and zeroth order for an electroweak decay signal. 

Our approximation for the jet cross sections is based 
on factorization for the relatively-collinear partons that 
form a jet from the remainder of process Q. For a jet of 
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cone size R, contributions that do not vanish as a power 
of R are generated by a function that depends only on 
the flavor of the parent parton, its transverse momentum, 
and the factorization scale. Denoting an jet shape by e, 
we then have, 



da 



dm j de 



dp. 



da c (p T ) dJ c {e,mj,p T ,R) 



dpi 



de 



(1) 

where da/dpT includes the hard scattering and the par- 
ton distributions of the incoming hadrons, and where the 
jet function for partons c in the final state is defined for- 
mally as in Refs. [(I [l3j]. 

In Ref. Q, we have found that the distribution of QCD 
jet masses in the range of hundreds of GeV is fairly well 
described by the jet function in Eq. {T]) at order a s , 
based on two-body final states. It thus seems natural 
to anticipate that for QCD jets, energy flow inside the 
cone would produce a linear deposition in the detector 
0, H, [ll|, [H, [l3|. While this is certainly the case for 
an event consisting of two sub-jets, it is a simpler con- 
dition, and more easily quantified. Indeed, such a linear 
flow should also characterize jets from the two-body de- 
cay of a highly-boosted, massive particle. We will see 
below that relatively simple jet shapes can help distin- 
guish QCD jets from many top-decay jets that involve 
three-body decay. We will also see that jet shapes can 
help separate samples that contain both QCD jets and 
jets from two-body decays, such as those of the W, Z or 
h. We emphasize that a single event may be analyzed by 
a variety of jet shapes, so that the resolution associated 
with each one need not be dramatic, so long as they are 
effectively independent. 

Top decay and planar flow. The linear flow of QCD 
jets at leading order should be compared with a > 3-body 
decay where the energy deposition tends to be planar, 
covering a two-dimensional region of the detector. An 
IR-safe jet shape, which we denote as planar floW) & two- 
dimensional version of the U D parameter" [H, H E3, 
[l8l |. distinguishes planar from linear configurations. The 
utility of a closely-related observable was emphasized in 
Ref. d3. 

Planar flow is defined as follows. We first construct for 
a given jet a matrix I w as 

1 \ - Pi.k Pi,l ,„s 

z^ Wi ~ — ' ( 2 ) 



jkl _ 



mj 



Wi Wi 



where mj is the jet mass, Wi is the energy of particle i 
in the jet, and pi^ is the k th component of its transverse 
momentum relative to the axis of the jet's momentum. 
Given I w , we define Pf for that jet as 



Pf 



4det(J ul ) _ 4AiA 2 
tr(J w ) 2 ~ (Ai + A 2 ) 2 



(3) 



where A1.2 are the eigenvalues of I w . Pf vanishes for lin- 
ear shapes and approaches unity for isotropic depositions 
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FIG. 1: The planar flow distribution for QCD and top jets 
obtained from MadGraph and Sherpa. Distributions are nor- 
malized to unit area. 



of energy. 

Jets with pure two-body kinematics have a differential 
jet function fixed at zero planar flow, 



I / dJ 



(4) 



This would apply to leading order for events with highly 
boosted weak gauge boson, Higgs and QCD jets. On the 
other hand, events that are characterized by > 3-body 
kinematics have a smooth distribution. 

Realistic QCD jets have, of course, nonzero Pf. Be- 
cause Pf is an IR safe observable, however, its average 
value can depend only on the hard momentum scales of 
the jet, that is, mj and px- This suggests an average 
Pf of order a s (mj) ~ 0.1 for high jet masses, times at 
most logarithms of that are order unity for these heavy 
jets. Correspondingly, higher orders corrections should, 
by analogy to two-jet event shapes (lbj |. replace the delta 
function of Eq. |(3} with a differential distribution that 
peaks near the origin and then falls off. For jets result- 
ing from three-body decay, on the other hand, we an- 
ticipate that corrections in a s shift the already-smooth 
distribution modestly, without affecting its overall shape. 
Finally, for the vast majority of high-p^ QCD jets, with 
masses mj <C Pt, planarity corrections associated with 
multi-gluon emission may be expected to be large, and 
to shift Pf to order unity. 

These considerations are confirmed in Fig. [TJ where 
we show the Pf distribution for QCD jet and hadronic 
tt events, for R = 0.4, p T = 1000 GeV and mj = 
140 - 210 GeV as obtained from MadGraph and 
Sherpa [2(| with jet reconstruction via (the IR-safe algo- 
rithm) SISCone [7|. We see that QCD jets peak around 
small values of Pf, while the top jet events are more 
dispersed. A planar flow cut around 0.5 would clearly re- 
move a considerably larger proportion of QCD jets than 
top jets. Correspondingly, we have confirmed by event 
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generator studies that low-mass QCD jets have much 
larger planar flows that those in Fig. [T] 

Two-body decay. While planar flow can help enrich 
samples with characteristically three- and higher-body 
kinematics, other jet shapes can also provide additional 
information on events with relatively low Pf. Here, we 
will still wish to distinguish the QCD background from 
highly boosted electroweak gauge bosons or Higgs [14{ as 
well as from top jets whose Pf happens to be relatively 
low. We begin with jets that are linear at lowest order, 
and identify a set of jet shapes that have some power to 
distinguish between the two. Fixing px, R and mj leaves 
only one free parameter to characterize the shape. 

The QCD jet function for two-body kinematics is de- 
fined as a matrix element in and is readily expressed as 
an integral over 9 S , the angle of the softer particle relative 
to the jet momentum axis. For a quark jet, for example, 
the integrand is therefore the differential jet function, 



a s C F f3 z z 2 



djQ CD _ 

d(c9 s ) ~ 4Trm 2 j(l - &c0,)(2(l - (3 z c9 s ) - z 2 ) 

(2(1 + &)(1 - /3 z c9 s ) - z 2 (l + c9 s )) 2 



z 2 (l + c6 s )(l - p z c6 s ) 
z\l + c9 s ) 2 



(1 - &c0 a )(2(l + p z )(l - P z c9 s ) - z 2 {\ + c9 s )) 



+ 3(l + /3 2 ) + 



(5) 



where 



= mj/px, (3 Z = \/l - z 2 and c9 s = 
cos 9 a . The jet mass function is obtained by the inte- 
gral d9 s (dJ/d6 s ), where 9 m is the angle with the 
smallest possible value of the softest particle, 9 m 
cos 1 (Vl — z 2 ), at which both particles have the same 
energy and angle to the axis. 

For signal events from a highly-boosted massive gauge 
bosons, we consider separately the cases when it is lon- 
gitudinal and when it has helicity (h = ±1), 



C 



dJ Lo„g ^ 

d(c9 s ) ~ (1 - x cB.) 2 
dJ h=±1 C 



(6) 



d(c9 s ) (1 - p z c9 s f 



1 - 



(zs9 s f 



2(1 -p z c6 s f 



where s9 s = sin0 s and C is a proportionality coefficient, 
totally determined from the two-body decay kinematics. 
We can interpret the appropriately normalized differen- 
tial jet functions, P X (9 S ) = (dJ x /d9 s )/J x as the proba- 
bility to find the softer particle at an angle between 9 S 
and 9 S + 59 s . As the ratio z decreases, the decay prod- 
ucts become boosted and the cone shrinks. For QCD 
jets from light partons, however, this shrinkage is much 
less pronounced. Plots of these jet functions show that 
the gauge boson distributions of Eq. ([6]) fall off with 9 S 
faster than do QCD jets, Eq. ((5]). This observation sug- 
gests that the signal (vector boson-jet) and background 
(QCD jets) have different shapes for fixed pt, R and jet 
mass. This may be used to obtain an improved rejection 



power against background events. We now consider a 
class of j et shapes, angularities, originally introduced in 
Ref. |l3l |2~0 ] for two-jet events in e + e~ annihilation. A 
natural generalization of these jet shapes to single-cone 
jets of large mass mj is 



To,{R,Pt) 



1 

mj 



E 

i£jet 



uji sin 



7T0A 

2R J 



2R 



(7) 

with mj the jet mass. The arguments of the trigono- 
metric functions vary from zero to tt/2 as 9 increases 
from zero to R, that is, over the size of the cone. These 
weights revert to the an gula rities as defined in for lep- 
tonic annihilation in 13j, |2l| when R — tt/2, so that the 
cone is enlarged to a hemisphere and mj is replaced by 
the center-of-mass energy in a two-jet event. For mas- 
sive jets, the angularities are clearly non-zero at lowest 
order, in contrast to the lowest order planar flow, Eq. 
dH) . Then, precisely because their IR safety, higher-order 
corrections to the r a distributions should be moderate. 

As the parameter a varies, the angularities give more 
or less weight to particles at the edge of the cone com- 
pared to those near the center. From the differential jet 
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FIG. 2: The ratio between the signal and background proba- 
bilities to have jet angularity f_2, R T ~ 2 . 

distribution functions in Eqs. © and ^) and the defi- 
nition of f a we can obtain the expression for P x (f a ) [as 
before x= sig (signal) or QCD] the probability to find a 
jet with with an angularity value between f a and f a + Sf a 
at fixed pT,R,mj and a. Our focus is not on the form 
of the individual distributions but rather on the ratio of 
the signal to background 



R(f a ) = 



P s[g (r a ) 

p Q CD (fa) 



(8) 



In Fig. [2] we show R Ta for a = —2 and z = 0.05 , for the 
different vector boson polarizations. In Fig. [3] we show 
the corresponding angularity distributions at the event 
generator level, comparing the output of MadGraph for 
longitudinal Z boson production to QCD jets in the same 
mass window. The pattern suggested by the lowest-order 
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Angularity (t. 2 ) 



FIG. 3: The angularity distribution for QCD (red-dashed 
curve) and longitudinal Z (black-solid curve) jets obtained 
from MadGraph. Both distributions are normalized to unit 
area. 



FIG. 4: Angularity, r_5 as a function of the azimuthical angle 
of the W(qq) pair, c/> q , for a typical top jet event, compared 
to the typical case two-body kinematics. 



prediction of Fig. [2] is confirmed by the output of the 
event generator, with signal and data curves crossing in 
Fig. linear r_ 2 = 0.02, where R(f- 2 ) ~ 1. 

Linear three-body decay. The leading-order differ- 
ential top jet function can be obtained by considering its 
three-body decay kinematics. The analytic expression is 
similar to Eq. ([6|) for the two-body case, although a bit 
more elaborate. In the following we simply point out 
a few features that may help angularities to distinguish 
top jets from background, even when they have relatively 
linear flow. 

The lowest-order three-body distribution is fully char- 
acterized by three angles. The first, 6b, is the angle be- 
tween the b quark and the jet axis. The second, 8y/ q , is 
the angle of the quark (from W-decay) relative to the W. 
The third, <fi, is the angle of the same quark relative to 
the plane defined by the W and the b. For an on-shell 
W, the distributions peak around 9b = m (as in two- 
body kinematics) and 6wq = @m(w) the minimal angle 
relative to the W momenta in the W rest frame. Because 
it is massive, the W's decay products move in somewhat 
different directions, even in the boosted frame, and their 
relative orientation induces the 0-dependence. Clearly, 
planar flow has maxima for odd multiples of <j) = 7r/2, 
and vanishes at lowest order at multiples of it. To tag 
top events at zero planar flow, angularities can be of use. 
In Fig. [4] we plot f_5 as a function of the azimuthical 
angle of the W(qq) pair, <p, for a typical top jet event. 
We also show the corresponding value for the two-body 
case (clearly <f> independent). For illustration we choose 
the kincmatical configuration that maximizes the corre- 
sponding differential jet distributions. We notice that 
this top angularity has maxima with </> at zero and ir at 
values far above the most likely two-body configuration. 
The reason is simply that angularities with large nega- 
tive values of a tend to emphasize flow at the edge of the 
cone. Other values of a weight individual jets differently 



in general. We consider this simple plot, along with the 
forgoing examples from event generators, as strong evi- 
dence for the potential of jet shape analysis. 

In summary, planar flow, angularities, and jet shapes 
that are as yet to be invented, may afford a variety of 
tools with which to distinguish the quantum mechanical 
histories of jets, whether resulting from heavy particle 
decay, or strong interactions. 
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